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Spectropolarimetry of distant sources of electromagnetic radiation at wavelengths 
ranging from infrared to ultraviolet are used to constrain Lorentz violation. A 
bound of 3 X 10""^^ is placed on coefficionts for Lorentz violation. 



Lorentz symmetry is an important part of our current understanding of 
particle physics. A violation of this symmetry would be a signal of physics 
beyond the standard modelEI For example, Lorentz violation can arise in string 
theoryo A general Lorentz-violating standard-model extension has been con- 
structedEl It allows for the possibility that the remnants of Lorentz violation 
occurring at the Planck scale may lead to small violations at energies attainable 
today. A number of experiments have been performed to test the fermion sector 
of the theory. The CPT-odd coefficients of the photon sector, which have been 
constrained experimentally to a highj-degree of precision, are expected to be 
zero from theoretical considerationsS'cl However, the CPT-even coefficients of 
the photon sector have received little attention. It is the goal of this work to 
understand the effects of these coefficients and to place bounds on them using 
existing spectropolarimetric measurements of distant cosmological sources. 

A Lorentz-violating electrodynamics can be extracted from the standard- 
model extension. We neglect the CPT-odd coefficients for the reasons men- 
tioned above. The relevant Lagrangian is £ = - jF^i^F^^'' - jikp) ^x^^yF"^^ F'^" 
where F^^ is the usual field strength, F^i, = d^A^ — d,jA^. The second term 
is CPT even and Lorentz violating. The coefficient {kp)i.^\^u is dimensionless. 
It has the symmetries of the Riemann tensor and a zero double trace. This 
leaves 19 independent components. 

The equations of motion resulting from this Lagrangian are modified in- 
homogeneous Maxwell equations. With the aid of the usual homogeneous 
Maxwell equations, plane-wave solutions of the form F^v{x) = F^^{p)e~^P"^ 
can be foundp'u The modified dispersion relation to leading order in {kF)K\^iv 
is 

P°± = (l-fp±fT)|p1, (1) 

where p = -^k^" and ct^ ^ ^{kafi)^ - p^, with fc"'^ = {kpT^'^^Pt^Pu and 
= l\p\. At leading order in {kp)^\f^^, the corresponding solutions for the 
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electric field, E±, are orthogonal and each is perpendicular to its group velocity 
Vg± = Vp- pj.. This implies the unit vectors e± = E±/\E±\ form a basis for 
the electric field. The general solution is of the form E{x) = (i?_|_e+e~*^+* + 
E-£-e~^^-*)e^P'^ . The fact that the phase velocities Vp± = p^p/p"^ of the two 
modes differ implies as the light propagates the relative phase between modes 
changes. The change in relative phase is 

A0 = (p° - )t w 27rAt;pL/A w AiTaL/X, (2) 

where L is the distance the radiation traveled and A is its wavelength. This 
phase change results in a change in the polarization of the radiation. The 
L/A dependence suggests, for very distant sources producing light at short 
wavelengths, tiny differences in phase velocity may become detectable. 

It is this L/\ dependence that is exploited in this work in order to ob- 
tain a bound on (fcF)KA/ii/- Recent spectropolarimptry of distant galaxies at 
wavelengths ranging from infrared to ultraviolet Q^IIj has made it possible to 
achieve values of L/A greater than 10'^^. Measured polarization parameters are 
typically order 1. Therefore ,we expect an experimental sensitivity of 10~^^ or 
better to components of {kF)K\nu- 

The first step in our analysis is to choose a coordinate system in which to 
work. A natural choose is a celestial equatorial system with the 3-axis aligned 
along the celestial north pole at equinox 2000.0 at a declination 90°. The 1- 
and 2-axis are at declination 0° and right ascension 0° and 90°, respectively. 
The goal is to place bounds on components of {kp)i^\f^u in this frame. However, 
for a source at an arbitrary position on the sky, this is not the most convenient 
coordinate system. Polarization is given by the behavior of E in the plane 
perpendicular to the direction of propagation. Therefore, for each source, we 
define a 'primed' frame where p''^ — (1; 0,0,1) at leading order. Then the 
primed-frame basis vector gg points from the source towards the Earth. To 
match standard polarimetric conventions we choose e'l so that it points south. 
The idea is to do much of the analysis in the primed frame where things are 
simple and then use observer covariance to write {kpYi^Xf^i, in terms of {kp)i^\f^i,. 
The two frames are related by a rotation. 

In the primed frame p = i(fc'" + fc'22) and cr^ ^ (fc'i2)2 ^ i (^ii _ ^22^2^ 

The form of suggests defining an angle ^ such that fc'^^ = cr sin ^ and 
— fc'22) = CTCOs^. Note that while p and <t are frame independent, 
^ is not. Solving the modified Maxwell equations in this frame gives i± cx 
(sin^, ±1 — cos^, 0). This implies that the birefringent modes are linearly po- 
larized. From the solutions e± and Eq. it is evident that a and ^ are the 
relevant parameters for polarimetry of a particular source. More precisely. 



2 



crsin^ and crcosf represent the minimal linear combinations of {kp)i^\p^u ef- 
fecting the polarization of a given source. It can be shown that asin^ and 
(Tcos^, written in terms of the celestial equatorial {kp),iX^^, depend on the 
right ascension and declination of the source and 10 independent components 
of {kp)f^x^i,. It is these ten components that are bounded in this work. We 
denote these components as fc", a = 1, ID. A suitable choose for k'^ in terms 
of (fcF)KA^i/ can be found in the literatureE 

In general, plane waves are elliptically polarized. The ellipse can be char- 
acterized by two angles: ip, the angle between e[ and the major axis of the 
ellipse and y = ± arctan , which describes the shape of the ellipse 

^ major axis ' ^ ^ 

and the helicity. The change in relative phase, Eq. (||), results in a change in 
both these angles. Most published polarimetric data of astronomical sources 
do not include measurements of x- Therefore, our analysis focuses on finding 
an expression for the change in -0. It will not only depend on the wave- 
length A, and the distance to the source L, but also on the values of ip and x 
when the light is emitted. Our approach is to look for wavelength dependence 
in the observed polarization. This approach assumes that the polarization at 
the source is relatively constant over the wavelengths considered. 

We seek an expression for S^p = ip — -tpo, the difference between ijj at two 
wavelengths, A and Aq. We find 

, sin|cosCo -f cos^sinCoCOs(50- (?!)o) 

1)7/; =2 tan = , (3) 

cos ^ cos Co — sin ^ sin Co cos{S(j) — 0o) 

where S(f) — A'k<7{L/\ — X/Aq), ^ ~ ^ — 2^o and 0o = tan^-^(tan 2xo/ sin^), 
Co = cos^^(cos2xoCOS^)H The polarization at Aq is given by V'o and xo- Two 
of these parameters need to be fit to the data. This is equivalent to fitting the 
initial polarization. The third parameter can be fixed to a convenient value. 

Table 1 lists 16 sources with published polarimetric data with observed 
wavelengths ranging from 400 to 2200 nm. For each source, we choose i/jq as 
the mean polarization angle and use Eq. (^) to create a x^ distribution. Each 
distribution is a function of ?/;, ^, Aq, and xo- They are then minimized with 
respect to Aq, and xo- 

Figure 1 shows the minimized distribution for the source 3CR 68.1. The 
features of this contour are common to all sources in Table 1. The contour 
corresponds to a confidence level of about (100 — 10~^)%. We see from Fig. 1 
that the parameter space away from tj = 0, ^ = 0°, ±90° are eliminated by this 
source. These are only the regions where the theory predicts no change in -0. 
The regions near ^ = 0°, ±90° correspond to the radiation being in a specific 
combination of birefringent modes. For example, ^ = 0° occurs if the light is 
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emitted in only one mode. We assume the probability of this happening for 
all 16 sources is small. With this assumption, the can be used to place a 
conservative constraint on a. In Fig. 1, the bound is shown as a horizontal 
line. 



Source 



a 



IC 5063Q p 
3A 0557-3830 
IRAS 18325-5925 
IRAS 13580-1818 
3C 324 
3C 256 
3C 356 

F J084044.5-h363328|i2 
F J155633_8-h351758 
3CR68.1EI 
QSO J2359-124iy 
3C 234Ei 
4C 40.361^ 
4C 48.480 
lAU O211-1220L 
lAU 0828+1930 
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L (Gpc) 



0.04 
0.12 
0.07 
0.13 
1.69 
1.92 
1.62 
1.71 
1.82 
1.70 
1.48 
0.55 
2.02 
2.04 
2.04 
2.08 



1030i/A 



0.56 - 2.8 
2.2 - 8.4 
1.0 - 4.9 



1. 



9.1 



58 - 130 
70 - 140 
57 - 120 
62 - 120 
67- 110 

59 - 130 
87 - 90 
51 - 75 

73 - 160 
75 - 160 

74 - 160 
78 - 160 



logi 



-30.8 
-31.2 
-31.0 
-31.0 
-32.2 
-32.2 
-32.2 
-32.2 
-32.2 
-32.2 
-31.1 
-31.7 
-32.2 
-32.2 
-32.2 
-32.2 



Table 1. Source Data. 




Figure 1: Contours of for the source 3CR 68.1. 



The bounds for each source are listed in the last column of Table 1. To es- 
timate the constraint on we assume, for each source, the data are consistent 
with (7 = 0. The bounds can be thought of as an estimate of the error b4> in a 
null measurement. We construct a second distribution, -y^ = Yjj{ajY/{8(jj), 
where the sum ranges over the 16 sources. This is a quadratic form in the 
fc" coefficients. A constant value of correspond to a ten-dimensional ellip- 
soid in the fc" space. We place a bound on the magnitude = Vk°lc^ by 
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minimizing with respect to the other nine degrees of freedom. This yields 
a bound of |fc°| < 3 x 10~^^ at the 90% confidence leveL 

These are the first bounds on the coefficients {kF)KXfiv They are compa- 
rable to the best existing bounds in the fermion sector of the standard-model 
extension. An improvement in this bound can be expected if more measure- 
ments similar to those used here are made. Similar measurements of x could 
also be used to impipve the bound. In the future it may be possible to include 
X-ray polarimetry^ which may lead to an improvement of several orders of 
magnitude. 
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